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Abstract 

We consider the scalar wave equation in the Kerr geometry for Cauchy data which 
is smooth and compactly supported outside the event horizon. We derive an integral 
representation which expresses the solution as a superposition of solutions of the radial 
and angular ODEs which arise in the separation of variables. In particular, we prove 
completeness of the solutions of the separated ODEs. 

This integral representation is a suitable starting point for a detailed analysis of 
the long-time dynamics of scalar waves in the Kerr geometry. 
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1 Introduction 



In a recent paper [8], the long-term behavior of Dirac spinor fields in the Kerr-Newman 
geometry, which describes a charged rotating black hole in equilibrium, was investigated. 
It was shown that solutions of the Dirac equation for Cauchy data in 1? outside the event 
horizon and bounded near the event horizon, decay in L?° as t — > oo. In this paper, we 
turn our attention to the scalar wave equation in the Kerr geometry. Our main result is 
to derive an integral representation for the propagator, similar to the one obtained for the 
Dirac equation in [8]. In our next paper [9], we will use this integral representation to 
analyze the long-time dynamics and the decay of solutions in . 

The analysis of the wave equation is quite different from that for the Dirac equation. 
The main difficulty is that, in contrast to the Dirac equation, there is no conserved density 
for the scalar wave equation which is positive everywhere outside the event horizon. This 
is due to the fact that the charge density, which was positive for the Dirac equation, 
is not positive for the wave equation. The other conserved density, the energy density, 
is non-positive either: it is in general negative inside the ergosphere, a region outside 
the event horizon in which the Killing vector corresponding to time translations becomes 
space-like. For these reasons, it is not possible to introduce a positive scalar product which 
is conserved in time. In more technical terms, we are faced with the difficulty that it is 
impossible to represent the Hamiltonian (i.e. the operator generating time translations) 
selfadjoint operator on a Hilbert space. 

We remark that the existence of the ergosphere is a direct consequence of the fact 
that the Kerr black hole has angular momentum [4]. Thus the ergosphere vanishes in the 
spherically symmetric limit. This simplifies the analysis considerably. 

A number of important contributions have been made to the rigorous study of the 
scalar wave equation in black hole geometries. The current last word on the stability of 
spherical black holes under scalar wave perturbations is the paper by Kay and Wald [14|, 
who proved using energy estimates together with a reflection argument that all solutions 
of the wave equation in the Schwarzschild geometry are bounded in More recently, 
Klainerman, Machedon, and Stalker [15] proved decay in L^ c of spherically symmetric 
solutions. These papers use the spherical symmetry of the Schwarzschild metric in an 
essential way. Whiting [21 J proved the absence of exponentially growing modes for the 
Teukolsky equation with general spin s = 0, |, 1, . . . (the case s = gives the scalar wave 
equation). Whiting's approach is based on interesting differential and integral transforms, 
which for a fixed angular momentum mode and fixed energy, convert the reduced ODEs 
into an equation admitting a positive conserved energy. Beyer [2] studied the wave and 
Klein-Gordon equations in the Kerr metric, using an approach based on C° semigroup 
theory. He proved that for each angular momentum mode, the Cauchy problem is well- 
posed, and he also obtained a stability result for the Klein-Gordon equation, provided that 
the mass parameter in this equation is sufficiently large. Finally, Nicolas [18] constructs a 
global solution for a non-linear Klein-Gordon equation in Kerr. 

Since the Hamiltonian cannot be represented as a selfadjoint operator on a Hilbert 
space, we are forced to employ methods which are quite different from those which we 
used in [8]. More precisely, the conserved energy gives rise to an indefinite scalar product, 
with respect to which the Hamiltonian is selfadjoint. By considering the system in finite 
volume with Dirichlet boundary conditions, we can arrange that the scalar product is 
positive on the complement of a finite-dimensional subspace. This allows us to use the 
general theory of Pontrjagin spaces [U [16]. In particular, the Hamiltonian is essentially 
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selfadjoint, and has a spectral decomposition involving a finite set of complex spectral 
points, which appear in complex conjugate pairs, together with a discrete spectrum of real 
eigenvalues. We write the projectors onto the invariant subspaces as contour integrals of 
the resolvent. In order to obtain estimates for the resolvent, it is useful to consider the 
Hamiltonian as a non-selfadjoint operator on a Hilbert space. This procedure also works 
in the original infinite volume setting, and we derive operator estimates which compare 
the resolvent in finite volume to that in infinite volume. Using these estimates, we can 
represent the spectral projector corresponding to the non-real spectrum as integrals over 
contours which are not closed and lie inside a region of of the form | Im u\ < c ( 1 + 1 Re uj \ ) ~ 1 
around the real axis. At this point, we make use of the fact that the scalar wave equation 
in the Kerr geometry is separable into ordinary differential equations for the radial and 
angular parts [I]. For the angular equation, we rely on the results of [TO], where a spectral 
representation is obtained for the angular operator, and estimates for the eigenvalues and 
spectral projectors are derived. For the radial equation, we here derive rigorous estimates 
which are based on the semi-classical WKB approximation. Using these estimates, we 
can express the resolvent in terms of solutions of the ODEs. Using furthermore Whiting's 
result that the ODEs admit no normalizable solutions for complex ui, we can deform the 
contours onto the real line. This finally gives an integral representation for the propagator 
in terms of the solutions of the ODEs with uj real. 

To be more precise, recall that in Boyer-Lindquist coordinates (t, r, i?, (p) with r > 0, 
0< , $<7r,0<( / 9<27r, the Kerr metric takes the form [U [12] 

ds 2 = gjk dx J 'x k 

= ^ (eft - a sin 2 dip) 2 -u(j£ + d#A - ^ (a dt - (r 2 + a 2 ) dtp) 2 (1.1) 
with 

U(r,tf) = r 2 + a 2 cos 2 , A(r) = r 2 -2Mr + a 2 , 

where M and aM denote the mass and the angular momentum of the black hole, respec- 
tively. We shall restrict attention to the case M 2 > a 2 , because otherwise there is a naked 
singularity. In the non-extreme case M 2 > a 2 , the function A has two distinct zeros, 



r = M - \Jm 2 - a 2 and n = M + \J M 2 - a 2 , 

corresponding to the Cauchy and the event horizon, respectively. In the extreme case 
M 2 = a 2 , the Cauchy and event horizons coincide, 

r$ = r\ = M . 

We shall consider only the region r > r\ outside the event horizon, and thus A > 0. 

In order to determine the ergosphere, we consider the norm of the Killing vector £ = Jj, 

: , A - a 2 sin 2 •& r 2 - 2Mr + a 2 cos 2 . . 

OijCe = 9tt = jj = jj • (1.2) 

This shows that £ is space-like in the open region of space-time where 

r 2 - 2Mr + a 2 cos 2 < 0, (1.3) 

the so-called ergosphere. It is a bounded region of space outside the event horizon, and 
intersects the event horizon at the poles i? = 0, w. 
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The scalar wave equation in the Kerr geometry is 



1 d 



-5 <9x 



5 



(1.4) 



where 5 denotes the determinant of the metric gij. In Boyer-Lindquist coordinates this 
becomes 



d ^d 1 / 2 2\ ^ ^ 
dr dr A \ dt dip 



sm 17 



tocos'!? 



Scosi? 



1 f ■ 2 a 9 5 

■TH asm + — 
sm w \ at dip 



$ = 0. (1.5) 



In what follows, we denote the square bracket in this equation by □ (although strictly 
speaking, it is a scalar function times the wave operator in (II. 4h ). 

A key property of the wave equation in the Kerr metric is that can be separated into 
ordinary differential equations by making the usual multiplicative ansatz 



-iu}t—ikip 



R(r) 0(0), 



(1.6) 



where w is a quantum number which could be real or complex and which corresponds 
to the "energy", and k is an integer quantum number corresponding to the projection 
of angular momentum onto the axis of symmetry of the black hole. Substituting (II. 6p 
into (|1.5p . we see that 

□* = (K,, fe + A,,fc)$, (1.7) 

where TZ U ^ and Au,k are the radial and angular operators 



—A— 

dr dr 

d 

<9cosi9 



1 

A 



((r 2 + a 2 )u + akf 



sin 2 $ 



d 



+ 



1 



d cos # sin 2 $ 



(auo sin 2 $ + k) A 



(1.8) 
(1.9) 



We can therefore separate the variables r and to obtain for fixed uj and k the system of 
ODEs 

Kw, k Rx = -xr x , A w , k e x = xe x , (1.10) 

where the separation constant A is an eigenvalue of the angular operator Au,k an d can thus 
be regarded as an angular quantum number. In the spherically symmetric case (i.e. a = 0), 
A goes over to the usual eigenvalues A = 1(1 + 1) of total angular momentum. Since the 
/c-modes are obtained simply by expanding the (/^-dependence in a Fourier series, we can 
in what follows restrict attention to one fixed /c-mode and omit the index k. We point out 
that for the A-modes the situation is more difficult because A as well as the corresponding 
angular eigenfunction 0>(i?) will in general depend on u. As a consequence, it is at this 
point not clear how to decompose the initial data into A-modes. 

We now reformulate the wave equation in first-order Hamiltonian form. The resolvent 
of this Hamiltonian will be one of the main ingredients in the statement of our main 
theorem. Letting 

* - ( i ) • < L11 > 
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the wave equation (j 1 . 5 j) takes the form 

id t ^ = H^>, (1.12) 

where H is the Hamiltonian 

a (5 



H = ( ° \ ) ■ (1-13) 



Here a and /3 are the differential operators 



d r Ad r - d cos $ sin 2 tf<9 cos + ( ^ r^rr ) d 2 



a = a sm if 



a 2 1 



A sin 2 <& 



A 

p _ -2a(' r2+ /» 2 - a ^h^)''(^-l) i8 ^ 

We now state our main result. 
Theorem 1.1 Consider the Cauchy problem 

□ $ = 0, (*,i$$)(0,ir) = * (af) (1-14) 

/or initial data £ Co ^ 7 "!' 00 ) x •S' 2 ) 2 which is smooth and compactly supported outside 
the event horizon, in the slicing associated to the Boyer Lindquist coordinates. Then there 
exists a unique global solution ty(t) = ($(t),idt$>(t)) which can be represented as follows, 

-^E e ~^E h > (/ ~ /_) ^ e ~ iUJt (Qk,n(w) Soo(u) *§)(r,t?) . (1.15) 
Here the sums and integrals converge in L 2 0C . 

In the statement of this theorem we use the following notation. The function is the 
kih angular Fourier component of *$>q, i.e. 

1 f' 2n 

We consider u in the lower complex half plane {Imw < 0}, and C £ is a contour which 
joins the points uj = —oo with u> = oo and stays in an e-neighborhood of the real line. A 
typical example is 

C £ = {x-iee~ x2 : x£K). 

C £ is the complex conjugated contour. Thus the integrals in (|1 . 15[) can be thought of as 
a "contour integral around the real axis" (see Figure [I]) , in analogy to the well-known 
Cauchy integral formula for matrices 



c 



-iAt 



1 



2m 



— ie-^{A-u)-^ 



c 



where A is a finite-dimensional matrix and C a contour which encloses the whole spectrum 
of A. For given u and k, the wave operator is a sum of a radial operator 1Zu,k an d an 
angular operator Au,k- As shown in [ID], the angular operator has for ui near the real line a 
purely discrete spectrum consisting of eigenvalues (A n ) nS ]N (see Lemma l2~Tj) . The spectral 
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Figure 1: Integration contours 



projectors onto the corresponding eigenspaces, which are one-dimensional, are denoted 
by Qk,n{^)- Furthermore, we write the wave equation in the Hamiltonian form (jl.121 
I1.13P and let S'oo(w) = (H — be the resolvent (in a suitable Sobolev space). The 
operator product Qk,n £00 can be expressed in terms of solutions of the reduced ordinary 
differential equations (see Proposition 15.41 with g = s and s as in Lemma l5.ip . Relying 
on Whiting's mode stability [21], we shall see below that the integrand in (11.15P is well- 
defined and holomorphic in the lower half plane, and thus the value of the integrals is 
indeed independent of the choice of C e . If the integrand were continuous up to the real 
axis, we could in (|1.15p take the limit e \ to obtain an integral along the real line. 
However, we do not know whether the integrand in (|1.15p is continuous on the real axis. 
Thus the integral in (|1.15p can be regarded as an integral over the real axis, with an 
"ze-regularization procedure" of possible singularities (if for instance the integrand had a 
simple pole at u>o S M, this would give rise to a (^-contribution at u>o). 

We point out that the global existence and uniqueness of solutions of the Cauchy 
problem can be obtained more generally in globally hyperbolic space-times (see e.g. |17j). 
The main point of Theorem ll.il is that we give an explicit decomposition of the propagator 
as a superposition of solutions of the ODEs which arise in the separation of variables. In 
particular, Theorem 11.11 shows completeness in the sense that the solutions of the coupled 
ODEs for real oj form a basis of the solution space. The explicit form of fll . 15H is useful 
for the study of the dynamics, because the time-dependence of is related by a simple 
Fourier transform to the w-dependence of the integrand in (|1.15p . which can in turn be 
analyzed by getting suitable ODE estimates [9|. 

We finally remark that the case of the wave equation for a scalar field minimally 
coupled to an electromagnetic field, 

g> h {V j - ieAj)(V k - ieA k ) $ = 0, (1.16) 

could be treated by similar methods in the non-extreme Kerr-Newman geometry, where 
now the metric is given by (|l.ip with 

U(r,-d) = r 2 + a 2 cos 2 ■d , A(r) = r 2 - 2Mr + a 2 + q 2 , (1.17) 

and the electromagnetic potential is 

Aj dx j = (dt - a sin 2 1? dip) , (1.18) 

where q denotes the charge of the black hole, and the parameters M, a, q satisfy the 
inequality M 2 > a 2 + q 2 . 
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2 Preliminaries 



In this section we briefly recall the variational formulation of the wave equation and the 
separation of variables. Furthermore, we bring the equation into a first-order Hamilto- 
nian form. Finally, we introduce and discuss scalar products which are needed for the 
construction of the propagator. 

The wave equation (jl.5p is the Euler-Lagrange equation corresponding to the action 



/OO fOO f'l /"7T 

dt dr (i(costf) / d<p£($,V<f>) , (2.1) 
-oo Jr\ J-l JO 

where the Lagrangian C is given by 

C = -A\d r ^\ 2 + ^\((r 2 + a 2 )d t + ad^\ 2 

- sin 2 1 8 cos — U | (« sin 2 dd t + <9 V )$ | 2 . (2.2) 
sin v 

According to Noether's theorem, symmetries of the Lagrangian give rise to conserved 
quantities. The symmetry under local gauge transformations yields that the vector field 

J k = -Im(¥Vjt$) , (2.3) 

called the (electromagnetic) current, is divergence free, and integrating the normal compo- 
nent of this current over the hypersurface t = const yields the conserved charge Q. More 
precisely, 

Q[$] = dr j ^(costf) JJ ^Q, 
where Q is the charge density 

Q = — $ 

_ j> 2 + a 2 ) 2 -/._ aid v $\ 2 . j aS /. fli , id^ 



A \ r 2 + a 2 J \ a sin 

Moreover, since the Kerr metric is stationary, the Lagrangian is invariant under time 
translations. The corresponding conserved quantity is the energy E, 

E[<f>] = [ dr f d(costf) / ^-£, (2.4) 
Jri J-l Jo 27 T 

where £ is the energy density 

- dC ^ t -C = ( (r ' + Q2)2 - a 2 sin 2 #) \d t $\ 2 + A |a r d>| 2 



<9$t V A 

+ sin 2 tf|d cos1? 3>| 2 + f—L--^) |<^$| 2 . (2.5) 
\ sin it A y 

One sees that all the terms in the energy density are positive, except for the coefficient of 
\d v Q\ 2 , which is positive if and only if r 2 — 2Mr + a 2 cos 2 > 0, i.e. outside the ergosphere. 
As a consequence, E is in general not positive. 
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Our analysis is based on a few properties of the angular operator A w , which we now 
state. For real lo, the angular operator Au clearly is formally selfadjoint on L 2 (S 2 ). 
However, this is not sufficient for our purpose, because we need to consider the case that 
lo is complex. In this case, A u is a non-self adjoint operator. Nevertheless, we have the 
following spectral decomposition, which is proved in |10j . 

Lemma 2.1 (angular spectral decomposition) For any given c > 0, we define the 
open set U C C by the condition 

\Imu)\ < - £ r . (2.6) 

1 1 l + |Reo;| v ' 

Then there is an integer N and a family of operators Q n (<^) defined for n £ NU {0} and 
lo £ U with the following properties: 

(i) The Q n are holomorphic in to. 

(ii) Qo is a projector onto an N -dimensional invariant subspace of A u . For n > 0, 

the Q n are projectors onto one- dimensional eigenspaces of A^ with corresponding 
eigenvalues X n (oj). These eigenvalues satisfy a bound of the form 

\X n (oo)\ < C(n) (1 + M) (2.7) 

for suitable constants C(n). Furthermore, there is a parameter e > such that for 
all n S N and lo £ U , 

|A„H| > ne. (2.8) 

(iii) The Q n are complete, i.e. 



n=0 

with strong convergence of the series in L 2 (S 2 ). 
(iv) The Q n are uniformly bounded in L 2 (S 2 ), i.e. for all n £ No, 

||Q n || < Ci (2.9) 

with c\ independent of to and n. 

If c is sufficiently small, c < e, or the real part of to is sufficiently large, |Reu>| > C(c), 
one can choose N = 1, i.e. A u is diagonalizable with non- degenerate eigenvalues. 

The proof of this lemma is outlined as follows. If |Imu;| is sufficiently small, the imaginary 
part of the potential can be treated as a slightly non-selfadjoint perturbation (see |13} 
Chapter 5, § 4.5]), giving a spectral decomposition into one-dimensional eigenspaces. On 
the other hand, for any fixed lo £ C, an asymptotic analysis of the resolvent (A u — A) -1 for 
large A (see e.g. [61 Chapter 12]) yields a spectral decomposition into invariant subspaces 
which for large |A| are one-dimensional eigenspaces. Thus the difficult point in the above 
lemma is to show that N and the constant c\ can be chosen uniformly in lo £ U . To this 
end, one must show that for real lo, the eigenvalue gaps of the selfadjoint operator Au 
become large as |cj| — ► oo. These gap estimates are worked out in |10] by analyzing the 
solutions of the corresponding complex Riccati equation. 
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After separation (|1.6p . the reduced wave equation takes the form 



dr dr 
d 



sin t?- 



d 



+ 



1 



(aw sin 2 i? + k) 2 



0. 



(2.10) 



d cos $ d cos $ ' sin 2 ?? ' 
Under the separation, the above expressions for the charge and energy densities become 



Q 
£ 



l$l 



1$ 



2\2 



(r 2 + a : 
A 

, | ( r 2 + a 2 ) 2 
A 



Rew 



ak 



r 2 _|_ Q 2 



21,2 



a 2 *; 



(r 2 + a 2 ) 2 



a 2 sin 2 i? |^Re u + 
-a 2 sin 2 ?? ( U 2 - 



k 



a sin 2 
k 2 



a 2 sin 4 ?? 



COS I? 



$1 



(2-11) 



(2.12) 



It is a subtle point to find a scalar product <., .> which is well-suited to the analysis of 
the wave equation. It is desirable to choose the scalar product such that the Hamiltonian 
H is Hermitian (i.e. formally selfadjoint) with respect to it. Since H is the infinitesimal 
generator of time translations, H is Hermitian w.r. to <., .> if and only if the inner product 
<\l/,Vl/> is time independent for all solutions Vl/ = ($,i5j$) of the wave equation. This 
can for example be achieved by imposing that <\I / , 1 I'> should be equal to the energy E 
corresponding to ^f. This leads us to introduce a scalar product by polarizing the formula 
for the energy, (12.41 12. 5ft . We thus obtain the so-called energy scalar product 



/oo p 1 r 

dr J d(cos i?) | 



[r 2 + a 2 ) 2 
A 



a 2 sin 2 ?? d t <$>d t & 



-a d r $> d r ®' + sin 2 •& d cos a os fl $' + 



1 



sin 2 1? 



-)d^d^' 



(2.13) 



where again ® = (®,idt&) and = ($', idt&'). If is a solution of the reduced system 
of ODEs (fLTUI) . can be written as = (* W ,A, ^*w,a) wi th '^('M?) = #a(?~)S a (??). 
Integrating by parts and dropping the boundary terms (which is certainly admissible when 
we consider the system in finite volume or when Vl/ has compact support), we can substitute 
the radial and angular equations into (|2. 13j) to obtain 



'(r 2 + a 2 ) 



dr 



ri 
2\2 



d(cosi?) 



A 



a z sin 2 i? (id t + cJ)$ $ w \ + 2ak 



r 2 + a 2 
A 



1 



a;, A 



(2.14) 



In the special case Vl/ = Vl/', this reduces to 



<*w,A,*«,A> = 2w 

'( r 2 + q 2 ) 2 
A 



dr 



ri 



d(cos-d) |<3? a 



a 2 sin 2 7? I Re a; + afc 



r 2 + a 2 
A 



1 



(2.15) 



By construction, the Hamiltonian is Hermitian with respect to the energy scalar product. 
However, the energy scalar product is in general not positive definite. This is obvious 
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in (|2.13p because the factor (sin -2 '!? — a 2 /A) is negative inside the ergosphere. Likewise, 
the integrand in (|2.15p can be negative because the factor ak in the second term in the 
brackets can have any sign. 

Apart from the energy, also the charge Q gives rise to a conserved scalar product. It 
is a natural idea to try to obtain a positive scalar product by taking a suitable linear 
combination of these two scalar products. Unfortunately, comparing (|2. 12j) and (|2.1ip 
one sees that it is impossible to form a non-trivial linear combination of Q and £ which 
is manifestly positive everywhere. One might argue that a suitable linear combination 
might nevertheless be positive because the positive term A|<9 r <I>| 2 + sin 2 {}\d cos ^\ 2 might 
compensate the negative terms. However, comparing fl2. 15f) with (12. lip , one sees that 
there is a simple relation between the energy scalar product and the charge, 

<*W,A,*a,,A> = 2wQ[tf WiA ] , 

making it again impossible to form a linear combination such that the integrand of the 
corresponding scalar product is everywhere positive. Stephen Anco showed that it is 
indeed impossible to introduce a conserved density for the wave equation which gives rise 
to a positive definite scalar product pp. We conclude that if we want to consider H as a 
selfadjoint operator, the underlying scalar product will necessarily be indefinite. 

But we can clearly consider H as a non-self adjoint operator on a Hilbert space, and 
this point of view will indeed be useful for the estimates of Section UJ Our method 
for constructing a positive scalar product is to simply replace the negative term —a 2 / A 
in (|2,13p by a positive term. More precisely, we introduce the scalar product (., .) by 

= J dr J d(cos^)| ^ r2 + ^ - a 2 w?& S jdj$d t & 

+Ad^d r & + sin 2 $d^$d cosi} & + _L-^*0 *' + (r ' + 0,2)2 $ • (2.16) 

sin v ' A J 

We denote the corresponding Hilbert space by Ti and the norm by ||.||. This norm dom- 
inates the energy scalar product in the sense that there is a constant c\ > depending 
only on the geometry such that the "Schwarz-type" inequality 

| | < ci H^ll H^'ll (2.17) 

holds for all € H. 

We finally bring the Hamiltonian and the above inner products into a more convenient 
form. First, we introduce the Regge- Wheeler variable u by 

du r 2 + a 2 d r 2 + a 2 d 

7k = A ' ]fr = A 5u ' ( ' ' 

The variable u ranges over (—00,00) as r ranges over (n,oo). Furthermore, we introduce 
the functions 

p = r 2 + a 2 - a 2 sin 2 •& _ A = (2.19) 

r l + a z 

» ^ --(^-^) <"<» 



p \ r z + a z 

2 1 „2 



* ~(^ + a 2 + -^- 2 ) ,2,21. 
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as well as the operator 



A 



<9u 9« 



A . 

~2~, 2 A 5 2 



2j,2 



a 2 fc 



r 2 _|_ a 2 



(2.22) 



where A52 denotes the Laplacian on the 2-sphere (recall that the parameter k is fixed 
throughout). Then, after integrating by parts, our inner products can on C 2 
written as 



<^l,^2> 
(*1,*2) 



-00 
00 



du / dcos t? p (^1 



du / dcos $ p 



A 
1 

A + 5 
1 



and the Hamiltonian takes the form 



H 



1 

,4 /? 



The functions p, (3, and 5 satisfy for a suitable constant c > the bounds 



1 



< 



P 



r 2 + a 2 



< c 



|*| < c. 



We abbreviate the integration measure in (|2.23p and (|2.24p by 

dp := p du d cos 1? . 



x S 2 ) 2 be 

(2.23) 
(2.24) 

(2.25) 



(2.26) 



3 Spectral Properties of the Hamiltonian in a Finite Box 

We saw in the preceding section that the energy inner product (|2,13p . with respect to 
which the Hamiltonian is formally selfadjoint, is in general indefinite. This fact remains 
true even when as in [8] we consider the system in a "finite box," i.e. when the range of the 
radial variable r is restricted to a bounded interval r G [r£,rjj] with n < Vl < tr < 00. 
Accordingly, in order to derive a spectral representation for the propagator corresponding 
to the wave equation (jl.5p . we will need to consider the spectral theory of operators on 
indefinite inner product spaces. Since there is an extensive literature on this topic, we 
here only recall the basic facts needed for our analysis, referring the reader to [3j [16] for 
details. 

A Krein space is a complex vector space tC endowed with a non-degenerate inner 
product <. , .> and an orthogonal direct sum decomposition 

/c = /c+e/c_, (3.1) 

such that (/C+, <., .>) and (ZC_, — <., •>) are both Hilbert spaces. A selfadjoint operator A 
on a Krein space K, is said to be definitizable if there exists a non-constant real polynomial 
p of degree k such that 

<p(A)x, x> > (3.2) 

for all x 6 T>(A k ). Definitizable operators have a spectral decomposition, which is similar 
to the spectral theorem in Hilbert spaces, except that there is in general an additional finite 
point spectrum in the complex plane (see [31 p. 180], [16, Thm 3.2, p. 34] and Lemma [3731 
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below). An important special case of a Krein space is when KL is positive except on a 
finite-dimensional subspace, i.e. 



K := dim/C_ < oo. (3-3) 

In this case the Krein space is called a Pontrjagin space of index k. Classical results of 
Pontrjagin (see [31 Thms 7.2 and 7.3, p. 200] and |16t p. 11-12]) yield that any selfadjoint 
operator A on a Pontrjagin space is definitizable, and that it has a K-dimensional negative 
subspace which is ^-invariant. 

We now explain how the abstract theory applies to the wave equation in the Kerr ge- 
ometry. In order to have a spectral theorem, the Hamiltonian must be definitizable. There 
is no reason why H should be definitizable on the whole space (r\, oo) x S 2 , and this leads 
us to consider the wave equation in "finite volume" [ri,rji] x S 2 with Dirichlet boundary 
conditions. Thus setting = ($, i<3?j) and regarding the two components (^i, ^2) of as 
independent functions, we consider the vector space V rL ,r R = (H ' 2 © L 2 )([ri,r^] x S 2 ) 
with Dirichlet boundary conditions 

*i(ri) = = fi(r fl ) . (3.4) 

Our definition of H 1,2 {[rL, tr] x S 2 ) coincides with that of the space W l,2 {{ri, tr) x 
S 2 ) in Section 7.5]. Note that we only impose boundary conditions on the first 
component ^1 of ^, which lies in H 1,2 . According to the trace theorem [7J Part II, 
Section 5.5, Theorem 1], the boundary values of a function in H ,2 {\tl, vr] x S 2 ) are in 
L 2 (S 2 ), and therefore we can impose Dirichlet boundary conditions. We endow this vector 
space with the inner product associated to the energy; i.e. in analogy to (|2.13p . 

rr R rl r/( r 2, 2)2 \ 

<^r, = J dr J d(cos &) I ^ — ^-i- - a 2 sin 2 tfj ^ 2 % 

+A^Td r *; + sin 2 $ d cos ^ 1 d cos ^[ + (-^ - ^) d^ld^'-X . (3.5) 

Lemma 3.1 For every r^ > r\ there is a countable set E C [r\,rji) such that for all 
tl £ (fi,rR) \ E, the inner product space V rL ,r R is a Pontrjagin space. The topology 
on V rLt r R is the same as that on (H ' 2 © L 2 )([r^,r^] x5 2 ). 

Proof. Since (|3.5|) involves no terms which mix the first component of with the second 
component, V rL ,r R clearly has an orthogonal direct sum decomposition V rL ,r R = Vi © V% 
with V1/2 = £ Pr L ,r R '■ *2/L = 0}- Furthermore, it is obvious that the space (V2, <., •> ) 
has a positive scalar product and that the corresponding norm is equivalent to the L 2 - 
norm. Hence it remains to consider V\, i.e. the space if 1,2 ([r£, Vr] x S 2 ) with Dirichlet 
boundary conditions and the inner product 

r-r R rl 

= / dr (i(cosi?) 

Jr L J-l 

x I Ad^d r & + sin 2 dcos*$' + (-^ - ^-J k 2 j .(3.6) 

Transforming to the variable u, (|2.18p . and using the representation (|2.23p . one sees that 
on the subspace C 2 ([ul,ur] x S 2 ) the inner product (|3.6f) can be written as 

<$ , = (<f>,A<f>') L 2 {[uL>UR]xS 2 >dlx) (3.7) 
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with A according to (|2.22|) . Here we set ul = u(tl), ur = u(rn), and dfi is the mea- 
sure (|2.26|) . A is a Schrodinger operator with smooth potential on a compact domain. 
Standard elliptic results [20, Proposition 2.1 and the remark before Proposition 2.7] yield 
that H is essentially selfadjoint in the Hilbert space TL = L 2 ([ul,ur] x S 2 ,dn). It has 
a purely discrete spectrum which is bounded from below and has no limit points. The 
corresponding eigenspaces are finite-dimensional, and the eigenfunctions are smooth. 

Let us analyze the kernel of A. Separating and using that the Laplacian on S 2 has 
eigenvalues —1(1 + 1), I £ No, A has a non-trivial kernel if and only if for some I G No, the 
solution of the ODE 



9,2 2^ 9 A a 2 k 2 



du v ; <9u ' r 2 + a 2 r 2 + a 2 



(u) = (3.8) 



with boundary conditions 4>(ur) = and 4>'(ur) = 1 vanishes at u = ul- Since this (j) has 
at most a countable number of zeros on (— oo,ur] (note that <f>(u) = implies (f>'(u) / 
because otherwise (ft would be trivial), <fi vanishes at ul only if ul £ Ei with Ei countable. 
We conclude that there is a countable set E = L)[Ei such that the kernel of A is trivial 
unless ul £ E. 

Assume that ul £ E. Then A has no kernel, and so we can decompose TC into the 
positive and negative spectral subspaces, H = TL+QTL-. Clearly, TC- is finite-dimensional. 
Since its vectors are smooth functions, we can consider TL^ as a subspace of V rL ,r R , 
and according to (|3.7j) it is a negative subspace. Its orthogonal complement in Vr L ,r R 
is contained in TL + and is therefore positive. We conclude that Vr L ,r R is positive except 
on a finite-dimensional subspace. 

It remains to show that the topology induced by < . , . > is equivalent to the H 1,2 - 
topology. Since on finite-dimensional spaces all norms are equivalent, it suffices to consider 
for any Ao > the spectral subspace for A > Ao, denoted by TC\ . We choose Ao such that 

a 2 k 2 

1 - A < \u ■ - min | — — — - } < 
Then for every $EC 2 n TC\ , 



[rt,r R ] V r 2 + a 2 



2 



(*) 

<*,^> = (y,A$>} L 2 {dfl) < c\\y\\ 2 Hl 

> I(*,^) L2(dM) + ^||*||| 2w 

> 1 m\\ 2 i v ° - 1 ihi'ii 2 4. A ° wmw 2 > 1 iim/ii 2 

where in (*) we used that the coefficients of the ODE ()3.8f) are bounded from above and 
below and that the zero order term is bounded from below by Vq. ■ 

We always choose tl and such that V rL ,r R is a Pontrjagin space and that our initial 
data is supported in [77,, r#] x S 2 . 

We now consider the Hamiltonian (II. 13ft on the Pontrjagin space V rL) r R with domain 
i^r] x S 2 ) 2 C 'Pr L) r R - For clarity, we shall often denote this operator by H rLiTR . 

Lemma 3.2 H rLirR has a selfadjoint extension inV Th . TR . 
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Proof. On the domain of H, the scalar product can be written in analogy to (|2.23p as 

= {^,S^) L 2 {[u ^ Ur]xS 2 M , 

where the operator S acts on the two components of ^ as the matrix 

where A is again given by (|2.22p and d/x is the measure (|2.26p . As shown in Lemma 13. 11 
S has a selfadjoint extension and is invertible. We introduce on C^((ul,ur) x S* 2 ) 2 the 
operator B by B = \S\~^ SH\S\~^ . The fact that H is symmetric in V TL) r R implies that B 
is symmetric in L 2 ([ul,ur] x S 2 ,dfi). A short calculation shows that 




\A\ \A\-?Ap \ 
p\A\-?A \A\+P 2 ) ' 



Treating the terms involving /3 as a relatively compact perturbation, we readily find 
that B 2 is selfadjoint on L 2 ([u L ,u R ] x S 2 ,dn) with domain V(B 2 ) = V(A)®V(A). Con- 
sequently, the spectral calculus gives us a selfadjoint extension of B with domain T>(B) = 
D(As) © V(A^). We extend H to the domain V{H) := \S\~* V{B). 

We now show that with this new domain, H is selfadjoint on V rL)TR . Suppose that for 
some vectors \£, G T > r L ,r R i 

i?V> = '> for all *' G X>(iT) . 

It then follows that 

{\S\H,B\S\H') L 2 {dll) = {S\S\-H,\S\H') L 2 W) for all *' G 

(note that the vectors G 7^,^ lie in H 1,2 © L? and thus both IS]^ S'|5'|~3$ e 
L 2 (dfi)). By definition of the domain of -ff, this is equivalent to 

(\S\^,B^) L2(dlM) = (S\S\-^,$) L2(dll) for all § G . 

Since -B is selfadjoint, it follows that the vector |S'|2\E' lies in the domain of B and that 
B \S\?V = S\S\~?$. This implies that * G £>(if) and that = * 



We now prove a basic lemma on the structure of the spectrum of the Hamiltonian 

Lemma 3.3 The spectrum of H rLyTR is purely discrete. It consists of finitely many com- 
plex spectral points appearing as complex conjugate pairs, and of an infinite sequence of 
real eigenvalues with no accumulation points. 

Proof. Since the operator H rL ^ TR is essentially selfadjoint, there exists a negative definite 
subspace L_ of T J rL ,r R of dimension k which is if rL ^-invariant (see [TU p. 11]). Let po 
denote the minimal polynomial of H rLTR on L_, i.e. 

Po{H rL ,r R )L- = 



14 



with degpo — K minimal. Furthermore, we let p be the real polynomial of degree < 2k 
defined by p = poPo- We claim that imp (H rL>rR ) is a positive semi-definite subspace. 
Indeed, we have for all x G Vr L ^ TR , 

<Po(H rLj r R )x,L-> = <x,p {H rL]TR )L_> = 0, (3.10) 

so that 

im P (H rLirR ) C imp {H rL>rR ) C (L_) x C (P rit ,r a )+, (3-11) 

as claimed. 

Next, since the square of the operator H rL>rR is elliptic, it follows that 

dimker (p (H rLirR ) < oo, (3-12) 

and from |16l Prp. 2.1], we know that for each eigenvalue £ of H rLi1 . R , the corresponding 
Jordan chain has finite length bounded by 2k + 1. It follows that p 2K+1 (H rLtJ . R ) has a 
finite-dimensional kernel and no Jordan chains. This implies that 

imp 2K+1 (H rLirR ) n ker/ re+1 {H rL)TR ) = {0} . (3.13) 

Furthermore, since the operator p 2k+1 (H rL ^ rR ) is selfadjoint, its image and kernel are 
clearly orthogonal. 

The image of p 2k+1 {H rL ^ R ) is contained in imp (H rL:1 . R ) and is therefore positive semi- 
definite. We shall now show that the space imp 2K+1 (H rLirR ) is actually positive definite. 
To this end, we let N be its null space, 

N :={x£ imp 2K+1 (H rL;rR ), <x , x>= 0}. (3.14) 

For all x G N and y G T>(H rLi r R ), we have 

<X, p 2K+1 (H rL ,r R )y>=0, (3.15) 

which is equivalent to 

<p 2K+1 (H rL , rR )x, y>=0, (3.16) 
because p is real. Since the scalar product is non-degenerate, this implies that 

p 2K+1 (Hr L ,r R )x = 0. (3.17) 

But we have just shown that kerp 2ft+1 (H rL rR ) and imp 2K+1 (H rLt1 . R ) have trivial inter- 
section. It follows that x = and therefore that imp 2K+1 {H rL)TR ) is positive definite, as 
claimed. 

Restricting H rL ^ R to im.p 2K+l (H rL! r R ), we have a selfadjoint operator on a Hilbert 
space. Thus the spectral theorem in Hilbert space applies, and the ellipticity of H 2 l Tr 
yields that the spectrum is purely discrete. On the finite-dimensional orthogonal comple- 
ment ker (p 2K+1 (H Ua>Ub )) we bring H rL ^ rR into the Jordan canonical form. ■ 
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4 Resolvent Estimates 



In this section we consider the Hamiltonian H as a non-selfadjoint operator on the Hilbert 
space Ti with the scalar product (., .) according to (|2.16p . We work either in infinite volume 
with domain of definition T>(H) = Cg°((ri,oo) x S 2 ) 2 or in the finite box r G [rL,rft] 
with domain of definition given by the functions in C°°((r^,r/j) x S 2 ) 2 which satisfy the 
boundary conditions (|3.4p . Some estimates will hold in the same way in finite and infinite 
volume. Whenever this is not the case, we distinguish between finite and infinite volume 



with the subscripts 



rL,rn 



and oo, respectively. We always consider a fixed A;- mode. 



The next lemma shows that the operator H — u> is invertible if either |Imu;| is large 
or |Imu;| 7^ and |Reu;| is large. The second case is more subtle, and we prove it using 
a spectral decomposition of the elliptic operator A which generates the energy scalar 
product. This lemma will be very useful in Section [71 because it will make it possible to 
move the contour integrals so close to the real axis that the angular estimates of Lemma [2.1l 
apply. By a slight abuse of notation we use the same notation for H and its closed 
extension. 



Lemma 4.1 There are constants c,K>0 such that for all E T)(H) and weC, 



(H - oj)^\\ > - I |Iuw 



A" 



1 + \Reuj\ 



1^1 



Proof. For every unit vector * £ T>(H), 



(H-uj)^\\ > !(¥,(# > |Im(*,(iJ-o>)*)| 
> |Imw| - -#*)*)!. 



(4.1) 



It is useful to work again in the variable u and the representation (|2.24p of the scalar 
product (., .) on C 2 (R x S 2 ) 2 . We introduce on Cq(R x S 2 ) 2 the operator H + by 



1 

A + 6 (3 



Comparing with (|2.24p one sees that H + is formally selfadjoint w.r. to the scalar product 
(.,.). Furthermore, one sees from (|2.25j) that H + differs from H only by a bounded 
operator, 



-5 



\H-H A 



< c. 



Thus on C 2 (R x S 2 ) 2 , 
\\H- H*\\ 



= \\(H - H+) - (H - H + 
< \\H-H+ II + \\(H-H 4 



2 WH-H, 



< 2c 



and substituting this bound into (|4.1|) . we conclude that 

\\(H > (\lmuj\ -c) 11*11 . 

In view of this inequality it remains to consider the case where |Reo;| is large. 

Using standard elliptic theory (see again [201 P- 86, Proposition 2.7] and [5])), the 
operator A with domain T>(A) = C^°(1R x S 2 ) is essentially selfadjoint on the Hilbert 
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space L 2 (dfi) := L 2 (IR x S 2 ,dfi), with dfi according to (|2.26|) . Clearly, A is bounded 
from below, A > — c, and thus a{A) C [— c, oo). For given A > 1 we let Pq and Pa be 
the spectral projectors corresponding to the sets [— c, A 2 ) and [A 2 ,oo), respectively. We 
decompose a vector W G W in the form = + with 

This decomposition is orthogonal w.r. to the energy scalar product, 

<#A,*o> = (*A, ( o 1 ) = °- 

However, our decomposition is not orthogonal w.r. to the scalar product (.,.), because 

(*A,*o) = (*A, ( M *q)v>(<M = (*A, ( o ) • 

But at least we obtain the following inequality, 

|(*A,*o)| < c||* || \\*l\\L*(dri , (4-2) 
where ^ A denotes the first component of ^a- Using that 

II*aI& ( *o = <*i,A- 1 *i>< ^I^aII 2 , 

we can also write (|4.2p in the more convenient form 

|(*A,*o)| < J ||*o|| II^All • 
Choosing A sufficiently large, we obtain 

||*|| 2 = I^aII 2 + 2Re(* A ,*o) + ||*o|| 2 < 4 (||* A || + ||*o||) 2 

and thus 

||*|| < 2(||tf A || + ||tf ||). (4.3) 
Furthermore, we can arrange by choosing A sufficiently large that 



2 



<*A,*A>= (*A, ( o J W)^) > \(**>yo 6 j)*A>xa(«^) = ~ ||*a| 

Next we estimate the inner products <*a,#*0>, (*o,^*a) and (^ ,H^ ). The 
calculations 

<¥ A ,ff*0> = (* A , ( 4 ) *0)i2 (d/l) = (* A , ( J? ) 

l(*o,^*o)| = (*o,(° ^)*o>l W 

< c ll*o||L2(d/i) + 2 11^*0 II L 2 (d/x) H^olU 2 ^) 
\\M\\h( d ,) = <*0, ( f ) *0>£ W 



< A 2 <* ,( J )*0>L^) = A 2 ||vl/ || 2 
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give us the bounds 

\<^A,H^ >\ < c ||^a|| ||*o|| 

|(^o, hv a )\ < cpoll ||*a|| 

\(* ,H9 )\ < (c + 2A) ll^oll 2 . 
Using the above inequalities, we can estimate the inner product <^A! (H — uj)^f> by 

| <^ A , [H - w)V> | > | <^a, [H - w)^a> | - | <*A, [H - uj)^o> | 

lima; I ,, , „o ,, T 
> J -^— 1 ||*a|| - c II*o|| ||^a|| • 

Applying the Cauchy-Schwarz inequality | <\I/ A , (H — o;) x I'> | < c\ H^aII ||(-ff — an d 
dividing by ||*a||) we obtain (possibly after increasing c) that 

> ^^II^aII - ||*o||- (4.4) 

c 

Next we estimate the inner product (\l/o, (H — a;)*), 

\(*o,(H-u)*)\ > \(* ,(H-u)*o)\-\(*o,(.H-u)* A )\ 

> (\u\ -c-2A) l^oll 2 " c ( X + x) A " ' 

We apply the Cauchy-Schwarz inequality (^q, [H — uj)*$>) < H^oll \\(H — u)fy\\ and divide 
by H^oll, 

\\(H - u)*\\ > (\u\ -c-2A) ||* || - c (\ + I^aII • (4.5) 

Choosing A = (|o;| — c)/4 and increasing c, the inequalities (|4.4j) and (|4.5|) give for suffi- 
ciently large \u\ the bounds 

\\(H-uj)^\\ > 1*^1 ||¥ A || - ||*o|| 

\\(H-u)*\\ > M||* ||- C ||* A ||. 

Multiplying the second inequality by 4/|o;| and adding the first inequality, we conclude 
that 

2[|(*- W )*[| > - ||* A || + ||*o|| . 

The result now follows from (14.31). ■ 



With the last lemma at hand, we are ready to introduce the resolvent. Namely, we let 



2K 

n = Iuj£C : \Jmu\ > : ; \ (4.6) 

1 + Reu; 



with K as in Lemma l4.1i 
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Corollary 4.2 If uj G fi, £/ie operator H — uj is invertible. The corresponding resolvent 

:= (H-uj)- 1 

satisfies the bound 



\\S{«>)\\ < (4.7) 
imw 



with c independent of uj G Q. 



Proof. In view of the preceding lemma, it suffices to show that the image of H — uj is 
dense in Ti for any ue!!. Otherwise, there would exist a non-zero £ Ti such that 

<(H - = for all $ G D(JT) , (4.8) 

that is a weak solution of the equation {H — uj)^ = 0. By the regularity theorem for 
elliptic operators on manifolds with boundary (cf. 19, Chapter5, Theorem 1.3]), every 
weak solution of this equation is a solution in the strong sense. These have been ruled out 
by the preceding Lemma. ■ 

Using (|4.6p in (|4.7p . we immediately get the bound 

< c(l + |Rew|) . (4.9) 

Since S(ui) is a bounded operator, its domain of definition can clearly be chosen to be the 
whole Hilbert space. We shall assume until the end of this section that uj G £1. 

The next lemma gives detailed estimates for the difference of the resolvents S rL)TR and 
Sqo in finite and infinite volume, respectively. By Q\(uj) we denote a given projector onto 
an invariant subspace of the angular operator A u corresponding to the spectral parameter 
A of dimension at most N (see Lemma 12.11 for details) . 

Lemma 4.3 For every G Co°((ri,r^) x S 2 ) 2 and every p G N, there is a constant 
C = C(^>,p) (independent of uj) such that 

|<¥, [S^,» - SooM] *>| < \^\ ■ (410) 

Furthermore, for every ^ G C^°((ri,r^) x S 2 ) 2 and every p G N and q > N, there is a 
constant C = C(^,p,q) (independent of uj and \) such that 

Proof. By definition of the resolvent, (H — uj) S(uj)^> = ^f. This relation holds both in 
finite and in infinite volume, and thus 

{{H - to) [S rL , rR (uj) - 5 M (w)] *) (r, 0) = if r L < r < r R . 

Iterating this identity and using the fact that H and S commute, we see that on [tl, tr] x 

" P+1 [Sr L ,r R (oj) ~ Soo(u)] * = [S rL ,r R (uj) ~ S^Uj)} H^ 1 * . (4.12) 
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Combining this identity with the Schwarz-type inequality (|2.17|) . we obtain 
|<tf, [S rL ,r R H ~ SocHl tf>| < ci \\S rL , rR (u) - ||*|| 2 

K+ 1 ||<^,[5 riirfl H- < s 00 H]^>| < I^^^H-^oo^)]^ 1 ^! 

< Ci \\S ri ,,r R (w) - S'ooMH || || . 

Since ^ is smooth and has compact support, H p+1 ^ also has these properties. The 
estimate (|4T9j) gives (|4.10p . 

In order to prove (|4.1ip . we first combine (|4.12j) with (|2.17p to obtain 

(1 + \to\ P+1 ) |<¥, Qx(S rL ,r R ~ ^oo)^>| 

< C! ||Q A || H^rK-^ooll 11*11 (11*11 + ||# p+1 \f||) ■ (4.13) 

Since g is at least as large as the dimension of the invariant subspace corresponding to A, 
(A, - A) 9 Q A = 0. Therefore, for every e C$°((r L ,r R ) x S 2 ) 2 , 

= <*, (A, - A) 9 Q A *'> = <{Al - A)«¥, Q A *'> • 

Expanding the power (.4* — A) 9 and using (|2.17p . we obtain 

iAr|<*,Q A *'>| < ^QiAr'H^^lliix^^QA*'!! 

1=1 

with combinatorial factors q (here X[r L ,r fl ] is the operator of multiplication by the char- 
acteristic function). Since the angular operator .A* is according to (jl.9p a polynomial in 
u; of degree two, the function (A^) l ^f is also polynomial in to, i.e. 

21 

(A*J* = £y* p , 

p=0 

where the functions are composed of ^ and its angular derivatives, as well as the 
coefficient functions of «4* . This gives the estimate 

21 

\\(A*jn < x>hi* p ii < c(i + n 21 ) 

with a constant c which depends only on \& and I. We thus obtain 

|A|*|<tf,Q A tt'>| < ^qW|A|^(1 + |o;| 2 0||xk^]Qa*'||. 
z=i 

Young's inequality allows us to compensate the lower powers of A, 

|A| 9 |<*,Q A *'>| < c(q,y)(l + \u\ 2 «)\\xir L ,r R ]Qxy'\\- 

We now choose ^' equal to the left side of (|4.12p with p = and p = r and take the sum 
of the resulting inequalities. Applying again the Schwarz inequality, we obtain 

|A|«(1+M r ) |<tt, Q A (S ri>rj? - ^oo)*>| < c(l+|u;| 2 «)||Q A || (||*|| + ||# r *||) . 

By choosing r sufficiently large, we can compensate the factor (1 + |oj | 2<? ) on the right. 
More precisely, 

|Ar(l + |^r 1 )|<M/,Q A (5 ri , rfl - < S 00 )M/>| < c / ||Q A ||||S rL , rii -5 00 ||(||^|| + || J ff p+2 '' +1 M/||). 
Adding this inequality to (|4.13p and substituting the estimate (14. 9p gives (14.111) . ■ 
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5 Separation of the Resolvent 



In this section we fix u cr(H), so that the resolvent S = (H — exists. As in the 
previous section, we assume that Q\ is a given projector onto a finite-dimensional invariant 
subspace of the angular operator Aw corresponding to the spectral parameter A. Our goal 
is to represent the operator product Q\S in terms of the solutions of the radial ODE. 
According to (fTTOj) and (TOD, the radial ODE is 



— A— 

dr dr 



(r 2 + a 
A" 



2^2 



0J + 



ak 



r 2 + a 2 



+ A 



R(r) 



(5.1) 



where A is the separation constant. We can assume that k > because otherwise we 
reverse the sign of u. We again work in the "tortoise variable" u, (|2,18p . and set 



0(r) = V ' r 2 + a 2 R{r) . 
Then equation (|5.ip can be written as 



1 d_ 

r 2 + a 2 du 

Using that 



(r 2 + a 2 



d_ 

du 



(r 2 + a 2 )^-(r 2 + a 2 ) 
ou 



+ \UJ + 



1 



ak 



r 2 + a 2 



AA 



(r 2 + a 2 ) 2 



\Jr 2 + a 2 



(5.2) 



(5.3) 



^{r 2 +a 2 ) 



(|5.3p simplifies to the Schrodinger-type equation 

d 2 



du 2 



+ V(u) <j>{u) 



(5.4) 



with the potential 
V(u) = 



ak 



r 2 + a 2 



+ 



A A 



(r 2 + a 2 ) 2 



+ 



Vr 2 + 



r l + a 1 



(5.5) 



We let 0i and 4>2 be two solutions of (|5.4h which are compatible with the boundary 
conditions. More precisely, in finite volume we satisfy the Dirichlet boundary conditions 
^1 (ul) = and 4>2(ur) = (again with ul = u(r^) and ur = u(rji)). Likewise, in infinite 
volume we only consider the case Im u < and let 0i and 4>2 be the fundamental solutions 
which decay exponentially at u = — oo and u = +oo, respectively (the existence of these 
fundamental solution will be established in Corollary 16. 4p . If the solutions 4>i and 02 were 
linearly dependent, they would give rise to a vector in the kernel of H — u>, in contradiction 
to our assumption uj &(H). Thus the Wronskian 



0i (u) 4> 2 (u) - 01 (it) 02 («) 



(5.6) 



is non-zero (note that w is by definition independent of u). 

We begin by constructing the "Green's function" corresponding to (15. 4)1 . 
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Lemma 5.1 The function 



s(u,u') := — TT~~I~\ >< I t U - U ', (5-7) 



satisfies the distributional equation 

d 2 



2 + V(u)\ s(u, v!) = 5(u - v!) . 
Proof. By definition of the distributional derivative, 

/•oo 

ri(u)(-dl + V)s(u,u')du = / ((-^ + V>(u)) s(u,u') du 



for every test function n £ Cq°(K). It is obvious from its definition that the function 
s(.,u') is smooth except at the point u = u' , where its first derivative has a discontinuity. 
Thus after splitting up the integral, we can integrate by parts twice to obtain 

((-d 2 u + V) V (u)) s(u,u') du 
u' 

rj(u) (—d 2 + V)s(u, u') du + lim n(u) d u s(u,u') 

oo 
/•oo 

+ / rj(u) (—d 2 . + V)s(u,u) du — lim rj(u) d u s(u, u') . 

J u i ' u\u' 

Since for u / u', s is a solution of (|5.4p . the obtained integrals vanish. Computing the 
limits with (|5.7p . we get 

((—d?, + V)r)(u))s(u,u)du = I lim — lim I n(u) d u s(u, u') 

' \u/u' u\u'J 

= i a 1 2 \ ^ u '^> ( < t > 'i( u ') <fo( u ') ~ ^K) M u ')) = , 

where in the last step we used the definition of the Wronskian (|5.6p . ■ 



In what follows we also regard s(u,u') as the integral kernel of a corresponding operator 
s, i.e. 



(s0)(n) := / du' s(u,u') 4>(u') du' . 

If Q\ projects onto an eigenspace of A^, we see from (|1.10p . (jl.7p . and (|5.2p that 

□ ((r 2 + a 2 yl Q x (i9, &) s(u, n')) = (r 2 + a 2 )-^ Q A (tf, #) 5(u - u') . (5.8) 

Loosely speaking, this relation means that the operator product Q\ s is an angular mode 
of the Green's function of the wave equation. Unfortunately, Q\ might project onto an 
invariant subspace of which is not an eigenspace. In this case, the angular operator 
has on the invariant subspace the "Jordan decomposition" 

AuQx = (A + AOQa (5.9) 

with J\f = J\f{uo, A) a nilpotent operator. Lemma 15.31 extends (|5.8p to this more general 
case. In preparation, we need to consider powers of the operator s. 
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Lemma 5.2 For every I £ No, £/ie operator s l is well-defined. Its kernel (s l )(u,u') has 
regularity C 2l ~ 2 . 

Proof. Writing out the operator products with the integral kernel, one sees that the 
operator s is obtained from s by iterated convolutions, 

s p+1 (u,u') = J s(u,u")s p (u",u')du" . (5.10) 

In the finite box, these convolution integrals are all finite because s(u,u') is continuous 
and the integration range is compact. In infinite volume, the function s(u, u') decays 
exponentially as u,u' — ► ±00 (see Corollary I6.4h . and so the integrals in (|5.10p are again 
finite. Hence s l is well-defined. 

Let us analyze the regularity of the integral kernel of s . By definition, s(u,u') is 
continuous, and (I5.10|) immediately shows that the same is true for s p (u,u'). Differen- 
tiating through (|5.10p and applying Lemma 15. 1\ one sees that s p satisfies for p > 1 the 
distributional equation 

This shows that incrementing p indeed increases the order of differentiability by two. ■ 



Lemma 5.3 For given A £ we let g be the operator 

00 

g = J2(-*0 l ( 5 - n ) 

1=0 

where N is the nilpotent matrix in the Jordan decomposition h5.9\) . Then 

□ ((r 2 + a 2 )-5 Q x (#, 0') g(u, u'fj = (r 2 + a 2 )^ Q A (0, 0') 5(u - u') . (5.12) 

Note that if Q\ projects onto an eigenspace, N vanishes and thus g = s. Furthermore, 
since M is nilpotent, the series in (|5.1ip is actually a finite sum. Thus in view of Lemma l5.2l 
(15. lip is indeed well-defined. 

Proof of Lemma \5.'A Denoting the radial operator with integral kernel 8(u — u') by l u , 
we can write the result of Lemma 15.11 in the compact form (— d 2 + V)s = l u . Hence on 
the invariant subspace, we can do a Neumann series calculation, 

00 00 
{-dl + V)g = Y^(-^) l (-d 2 u + V)s l+1 = J>A0^ = l u - M g , 

1=0 k=0 

to obtain that (-S 2 + V + Af) g = l u . According to (fTTOj) . (ITTfl) . and (^2|) . this is equiv- 
alent to (I5TT21 . ■ 



We come to the separation of the resolvent. In order to explain the difficulty, we point 
out that H and Q\ do not in general commute, and thus 

(H-u)Qx ^ Qx(H-u) and Q x S ^ S Q x . 
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Therefore, one must be very careful with the orders of multiplication; in particular, it 
is not possible to simplify the operator product (H — uj)Q\S. However, we know from 
the separation of variables that for every solution ^ of the equation (H — lv)^ = 0, its 
projection Qx(u))^f is again a solution. In other words, H and Q\ do commute on the 
kernel of {H — uj). This fact will be exploited in the proof of the following proposition. 

Proposition 5.4 For uj o~(H) we let Q\ be a spectral projector of the angular operator 
Aw- Then the resolvent of H has the representation 

Q x S(u) = QxT(uj,\) , 
where T is the operator with integral kernel 

T(u, 0; u', = £(cos 0- cos 0') <*(«-«') 

+5(cos0-cos0O (r 2 + a 2 ) *g(u,u')[ \ ,L \ , L • 5.13 

Here g is the Green's function h5.11\) (which depends onto and\), andr,a are the functions 

a = (r 2 + a 2 )-f ((r 2 + a 2 ) 2 - A a 2 sin 2 tf) 
t = 2ak (r 2 + a 2 )~i ((r 2 + a 2 ) - A) + uia . 

Proof. Let us compute the operator product (H — u>)QxT. We first consider the first sum- 
mand in (|5,13p . which we denote by T\. In this case, the operator product is particularly 
simple because the second column in the matrix (j!.13[) involves no u- or ^-derivatives. We 
obtain 

QaTOCu^u', 00 = Q A (0, #)6(u,u') (^l ) _ u o)- (5-14) 

Next we consider the second summand in (|5.13p . which we denote by T 2 . Fixing tt',0' 
and considering T 2 as a function of u, 0, we see from Lemma 15.31 that each column of 
Qa^2 is f° r u < u' & vector of the form = with $ a solution of the separated 

wave equation 111. 7ft . The same is true for u > u'. Hence for u ^ u 1 , Q\Ti is composed of 
eigenfunctions of the Hamiltonian, 

((F-w)QaT 2 )(u,0;«',0O = iiu^u'. 

It remains to compute the distributional contribution to (H — uj)QxT2 at u = u'. Since 
T 2 is continuous at u = u' , we only get a contribution when both radial derivatives act on 
the factor g. According to Lemma [5.21 the higher powers of s are in C 2 , and thus we may 
replace g by s. Applying (|1.13p . (|2.18p . and Lemma I5TTI we obtain 

)^T 2 )(.,0;<0O = ^Qa(0,0O^-O a(n % ). (5-15) 

We add (|5.14|) to (|5.15p and carry out the sum over A E a(Acj). Since the spectral 
projectors Qx are complete (see Lemma |2~T1 (iii) ) . YlxQxi 1 ^^') gives a contribution only 
for = 0'. We thus obtain a multiplication operator, 

(H-u>)Q x T(\) = l+^-t^ + I) 
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Using the explicit form of the functions r, a, and 0, one sees that the second term vanishes. 
Thus _ 

(H-u)Q x T(X) = 1. 

Multiplying from the left by QyS and using the orthogonality of the angular spectral 
projectors gives the result. ■ 



6 WKB Estimates 

In this section we shall derive estimates for the radial ODE (|5.4p in the regime 

|Reo>| 3> 1 and |Ima;| < c. (6-1) 

In this "high-energy regime", the semi-classical WKB-solution should be a good approx- 
imation. In order to quantify this statement rigorously, we shall make an ansatz for (j) 
which involves the WKB wave function and estimate the error. 
Our first lemma gives control of the sign of Re VV. 

Lemma 6.1 There is a constant C such that for all uj with 

Imw / , |Rew| > C 
and all A G a(A u) ), the function ReW has no zeros. 



Proof. At a zero of Re v V, the function V is real and non-positive. Thus it suffices to 
show that the imaginary part of V has no zeros. 

We first estimate the imaginary part of the angular spectrum. For any A G cr^A^) we 
let <3?a be a corresponding eigenvector. Then 



^(($,^) i2 -(M$) i2 ) = I($,(^-^)$) L2 

2i 2i 



Im(A) <<D,cD) L2 

where (., .)L 2 is the L 2 -scalar product on S 2 . Hence, according to (jl.9j) 



|ImA| < 



sup 

s 2 



Im 



1 



sin $ 



(ow sin 2 ?? + kf 



< 2a |Rea;| |Imu>| + |2a/clma;| 
The imaginary part of (15. 5p is computed to be 

ak 



ImV 



2 Reco + 



r 2 + a 2 



Imti) + 



A 



(r 2 + a 2 ) 2 



Im A 



(6.2) 



(6.3) 



Using (|6.2p . the second summand is estimated by 



A 



Im A 



< 2 



a 2 A 



(r 2 + a 



2\2 



Rewl H | |Imu;| 

a 



(r 2 + a 2 ) 2 

Z 2 ) 2 ViUUrvlU Mil MIA V VV LU 1 

smaller than one. Thus there is a constant c with a 2 A(r 2 + a 2 )~ 2 < c < 1. This shows 



The factor a 2 A(r 2 + a 2 ) 2 vanishes on the event horizon and at infinity and is always 
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that after choosing |Recj| sufficiently large, the first summand in 
second, and so ImV has no zeros. 



dominates the 



In what follows, we assume that the assumptions of the above lemma are satisfied. We 
choose the sign convention for the square root such that 



Re y/V(u), ReV(«)* > for all u £ 
Furthermore, we shall restrict attention to u in the range 

— c < Imw < 0, |Rea>| > C, 



(6.4) 



(6.5) 



where c is any fixed constant and C will be chosen depending on the particular application. 
The next lemma, which we will need in Section estimates yV inside the "finite box" 
[uli ur] uniformly in u for large Reuj. 

Lemma 6.2 For every angular momentum mode n and every c, e > 0, there are constants 
C and d such that for all u £ [ul,ur] and all uj in the range 16. 5\) . 



\ReW\ < c 
|Im y/V(u) - Im y/V(Beu)\ < e. 



(6.6) 
(6.7) 



Proof. We set ojq = Reu, A = A(wo) and introduce for a parameter r G [0, 1] the potential 
V T = V(u ) + tW with 



W 



-2i Ima; ( Rea; + 



a A: 



+ (Imw) 2 + (A-A ) 



A 



(r 2 + a 

Then Vq = V{ojq) and V\ = V(u). The mean value theorem yields that 

W 



2\2 



Re y/V{u) - Re y/V(u ) 
Im y/V(u) - Im ^V(lo ) 



< sup Re 

re[0.1] V 2 V^r 

T ( W 

< sup Im I — -j= 

re[0.1] VVK 



(6i 
(6.J 



By choosing C sufficiently large, we can clearly arrange that V(ujq) < 0, and thus 
Re ^JV(uiq) = 0. Furthermore, one sees immediately from the explicit formulas for V, 
W together with the estimate for the angular eigenvalue (|2.7|) that 



Re VF = O(|Rew| ) 
V^-iReu; = O(|Rew| ) 

Using this in (|6.8p and (|6.9[) gives the claim. 



ImW = OdRewl 1 ) 



We introduce the WKB solutions a and a by 



a( u ) = cV* exp W^j , a(u) = cFi exp ^-jT v 7 ! 7 ^) , 



(6.10) 
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where c and c are some normalization constants. A straightforward calculation shows that 
these functions satisfy the Schrodinger equation 

IV" ^ /v'\ 2 
a" = Va with V = V-- — + — { — \ . (6.11) 



We can hope that a and a are approximate solutions of the radial equation (|5.4p . In order 
to estimate the error, we first write fj5.4j) as a first order system, 

*'=(°J)* » ith •-($)■ <6 ' 12) 

Next we make for \E' the ansatz 

= A$ with A = ( a , ?,) (6.13) 

y a a J 

and $ a 2-component complex function. ^4 is the fundamental matrix of the ODE (|6.1ip 
and thus 

A'=(£J)A. (6.14) 
Differentiating through the ansatz for f)6. 13|) and using (|6.12l 16. 14|) , we obtain that 

A*'=( v °_ v (6.15) 

The determinant of A is a Wronskian and thus constant. A short computation using (|6.10p 
shows that det A = —26c. Hence we can easily compute the inverse of A by Cramer's rule, 

A -i = L ( ~" 

26c \ —a' a 

and multiplying (|6.15p by A' 1 gives 



1 /„ f>N ( —aa —a 2 s 
®=-W^( V - V )\ x2 ax 



26c V a aa 

Finally, we put in the explicit formulas (|6.1ip and (|6.10|) to obtain the equation 

$' = W ( ~ X 1 \ $ , (6.16) 



where W and / are the functions 

J_ (— _ - — 

86b \vJ ~ 4\4 



IT" - — ( — - - ^ , / = exp (2 j^Vv^j . (6.17) 

We shall now derive an estimate for the solutions of the ODE (|6.16p . The main 
difficulty is that when the function / is very large or close to zero, the matrix in (|6. 16j) 
has large norm, making it impossible to use simple Gronwall estimates. Instead, we can 
use that according to (16. 4p . the function 



I /| = exp (^2 £ Re W) 



is monotone. 
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Theorem 6.3 Assume that the potential V in the Schrddinger equation J5,-^[ j satisfies the 
conditions and that the function W defined by is in L X (R). Then there is a 

solution <I> of the system of ODEs \6. 16\) with boundary conditions 



lim = 

u — * — OO 

This solution satisfies for all u £l the bounds 



$i(u) — exp 



W 



< e 



4||VK|| 



\W\ 



|$ 2 («)| < e 4 ^ 1 H^llx |/(«)| 
Proof. We set /? = |/| and introduce new functions a and 6 by 

a = <&! and 

According to (|6.16|) . they satisfy the following ODEs, 



b = ^ 
P 



(6.18) 



-Wa-Wjb 
b' + ^-b = W-a + Wb. 



a 
j 



P P 

This gives rise to the following differential inequalities, 



—s/m = —Re (a a') < -\a\ReW + \b\ 
du \a 



W 



f 



\b\' 



,j 



^Re(6 6') < -\b\?- + 
\b\ P 



W 



+ \b\ ReW . 



Using that p is monotone and that |/| = p, we obtain the simple inequality 

(H + |6|)' < 2\W\(\a\ + \b\). 
Integrating this inequality from v to u, — oo < v < u < oo, gives the "Gronwall estimate" 

(|a| + \b\)(u) < (|a| + \b\)(v) exp ( 2 T |wf\ < (|a| + \b\)(v) e 2 W w ^ . (6.19) 



We now let be the solution of (|6.16p with boundary conditions $>( v \v) = (1,0). 
In order to estimate & v \ we rewrite (16.161) as 



f: w ^\u))' = We-K w f^K 
We integrate and use (|6.19|) to obtain the inequalities 



(«) 



< e 

< e 



3\\W\\ 



\w\ 



pW < e 



3IIW 



ru 

11 PW) / 



w\ 



(6.20) 
(6.21) 
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where in the last step we used the monotonicity of p. 

The inequalities (|6.20p and (|6.2ip yield that for every e > 0, there is a u such that for 
all v, v' < u, the exponential on the left side of (|6.20p and (|6.2ip are arbitrarily close to one, 
and the integrals on the right can be made arbitrarily small. Thus — $>( v >)(u)\ < e. 

Due to the factor p(u) on the right of (|6.2ip . we even know that |(ftC«) _&(«')) < e (with 
b according to (|6.18p ). Since (|6.16p is linear, — ) is also a solution. Applying (|6.19p 
for this solution and choosing v = u, we obtain that for all u > u, |(<&W — & v ^)(u)\ < ce 
with a constant c being independent of u. This shows that $( v '(u) converges as v — ► — oo, 
and that the above estimates are still true for v = — oo. 

The theorem now follows from (|6.20p and (|6.2ip if we set v = — oo and pull out a factor 
of e J-oo and e~ J -™ Respectively. ■ 



The above theorem has two immediate consequences: First, it yields the existence 
of solutions (ft and (ft which decay exponentially at minus and plus infinity, respectively. 
Second, it gives very good control of the global behavior of these solutions if |Recj| is 
large. 

Corollary 6.4 For every angular momentum mode n and every ui with Imcj < 0, there are 
solutions 4> o-nd 4> of the Schrddinger equation |5.^[ j which satisfy the boundary conditions 

lim 



= 1 = lim e iuju (j)(u) 



u— >oo 



Proof. It suffices to construct <p, because 4> is obtained in exactly the same way if one 
considers the ODEs backwards in u (i.e. after transforming the radial variable according 
to u — » —u). We choose $ as in Theorem 16.31 and let <^> = (j> be the corresponding solution 
of the Schrodinger equation given by (|6.13p and (|6.12p . Note that the corollary only makes 
a statement on the asymptotic behavior of <p as u — > — oo, and thus the behavior of (ft on 
any interval [«o,oo), < is irrelevant. Thus we may freely modify the potential V on 
any such interval. In particular, we can change the potential V on [uo, oo) such that it is 
constant for large u. For any e > 0, we choose uq so small and modify ^ U0iOO ) such that 
ReV^ > and ||W||i < e/3 (this is possible because V" decays for large \u\ at least at 
the rate ~ |ti|~ 3 ). Then Theorem 16.31 applies, and we obtain that 

|*i -1| < e, |$ 2 | < e\f\. 

Using these bounds in (16. 13f) . one sees that \<ft/a — 1| < e and thus, after choosing the 
normalization constants c and c in (|6.10p appropriately, 

limsup|e-^ u ^| < l + e and liminf \e~ iuJU <j>\ > 1-e. 

u — oo u >— 00 

Since e is arbitrary, the result follows. ■ 



Proposition 6.5 For every n and c, e > 0, there is a constant C > such that for all 
uj in the range $6. 5|) . the solutions (ft and (ft of Corollary \6.4\ are close to the (suitably 
normalized) WKB wave functions a and a, \6.10\) . in the sense that for all u E M, 













1 


< e and 




+ 








a 




a 
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The reason why we need to choose C large is that the functions y'/|y| 3//2 and W must 
be sufficiently small. More specifically, one can choose C such that 



\W h < 3 



and 



W'\ 



1 

V 2 3 



Proof of Proposition \6.5[ Using (|2.7p in (|5.5p . one sees that in the strip —2c < Imcj < 0, 
the potential V satisfies the bound 

\V(u)+uj 2 \ < d(l + |w|) . 



On the other hand, differentiating (15. 5h and using (12.18P and (12. 7j) . one sees that in this 
strip, 

+ < (i + M)^(«), 

where g is a function which decays for large |u| at least at the rate ~ u~ 2 . Putting these 
estimates for V and V" into (|6.17p . one sees that by choosing C sufficiently large, we can 
arrange that for all u in the range (|6.5p . ||W||i < e/3. Theorem 16.31 yields that 



1*1 -1| < e, |* 2 | < 
Dividing the first row in (|6.13p by a, we obtain the identity 



(6.22) 



a 



and using (|6.22j) gives 



a 



*1 + - *2 

a 



< e l + l/l 



^From the second row in (|6.13|) we obtain similarly, 



a 



1 



< e l + l/l 



a 



a' 



Finally, we apply the elementary estimates for the WKB wave functions 

2 



a 


1 




1 






a 




d' 


~ m 




5 1 
V~3 - Vi 



< 



l/l ' 



where in the last step we applied the above bounds for V and V and possibly increased 
C. 

The solution <fi is obtained similarly if one considers the Schrodinger equation (pT4 
backwards in u and repeats the above arguments. 



The next two propositions give estimates for composite expressions. 
Proposition 6.6 Under the assumptions of Proposition W. 

w(<j>,<f>) 



1 



< As. 
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Proof. Rewriting the Wronskian as 



6c' a — - — a a' 
a a 



a a 



we can put in the estimate of Proposition 16.51 to obtain 

w((j), cj>) — w(d, d) < 4e (\d! d\ + \d . 

Furthermore, a short explicit calculation using (|6.10p shows that 
w(d,d) = 2y/V act 
(da)' = -~6h = -\v'V-l W (d,d) 

and thus 

±(w(d,d) + (dd)') = ^(d,a) (l-^'Hj 



a a 



a a 



(w(d, d) — (d a)') 



Substituting these relations into (|6.23p gives 

w(<j),4>) 



w(d, d) 



< 4s[l + l 



(6.23) 

(6.24) 
(6.25) 



Here the left side only involves Wronskians and is thus independent of u. Hence we may 
on the right side take the limit u — > oo. This gives the result. ■ 



For ul < u R we set 



[U L ,U R \ 



a 



4>{u L ) 4>{u R ) - <j>(u R ) 4>(u L ) 

d(u L ) d(u R ) - d(u R ) d(u L ) 



Proposition 6.7 Under the assumptions of Proposition {6. 









< 8e exp (^2 j * Re VV^j 


sin ^2 J ImVfj 


a [u L ,u R ] 







-1 



Proof. Rewriting <f>[ UL)UR ] as 

?[u L ,u R ] = ot(u L ) a(u R ) — — r- — — r- - a[u R ) a[u L ) —, — r T7 — r 
L m a(u L ) a[u R ) a[u R ) a[u L ) 

Proposition 16.51 yields that 

\4>[u L ,u R ] - <*[u L ,u R ]\ < 4e(\d(u L ) d(u R )\ + \d(u R ) d(u L )\) 
Furthermore, it is obvious from (|6.10p that 



d(u L ) d(u R ) = d(u R ) d(u L ) exp 



■.rw) 

Jul / 



(6.26) 



(6.27) 
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and thus from (|6.26[) . 

a [u L ,u R ] = ot(u R ) a(u£) ^exp 

= a(u L ) a(u R ) ^1 - exp ^2 
Dividing (16.271) by (*[u L ,u R ] an d putting in the last identities, we obtain 

8e 



"L 



(6.28) 



[U L ,U R \ 



1 



a 



< 4e 



1 + 


exp 


r 2 C' 


/7) 




exp 


-i 



< 



where in the last step we used that Re v > 0. We finally estimate the obtained denom- 
inator from above, 



exp 



W)-i 



> 



Imexp [-2 1 W 



exp ^-2 ^ R Re v 7 ^ sin ^2 J " Im v 7 ^) 



7 Contour Deformations 

In this section we shall use contour integral methods to prove the main theorem. Recall 
that in Section [3l we showed that the Hamiltonian H ULtUR in finite volume is a selfadjoint 
operator on the Pontrjagin space V UL ,u R - It has a purely discrete spectrum, and for 
each to £ a(H ULtUR ), the projector E u onto the corresponding invariant subspace can be 
expressed as the contour integral 

K = -J- f Su L ,u R (u') du>' , 

2m J\u'-uj\<e 

where e is to be chosen so small that B £ {uj) contains no other points of the spectrum. 
The theory of Pontrjagin spaces also yields that a(H ULjUR ) will in general involve a finite 
number of non-real spectral points, which lie symmetrically around the real axis. We let 
Eq be the projector onto the invariant subspace corresponding to all non-real spectral 
points, 

E C := Yl E "- t 7 ' 1 ) 

LU£cr(H UL , UR )\R 

Our first lemma represents Eq as a Cauchy integral over an unbounded contour. More 
precisely, we choose a contour C UL)UR in the lower half plane which joins the points +oo 
with — oo and encloses the spectrum in the lower half plane from above. Furthermore, if 
Re a; is outside the finite interval [a>_,u;+], u should be in the open set Q (see (|4.6|) ) and 
should approach the real axis as |Imcj| ~ — |Retj| _1 (see Figure [2]). 
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Lemma 7.1 The spectral projector corresponding to the non-real spectrum (7.1) has the 
representation 



— lim Im 

7T L^oo 



L 3 



a (L + iuo) 3 
Proof. The Cauchy integral formula yields that 



(7.2) 



1 

2ttI 



S Ul ,u r {u) duj + 



c 



1 

2TTi 



S Ul , Ur (uj) (ko 



c 



7T 



Im (h S Ul , Ur {lu) doj , 



c 



where C is a closed contour which encloses the spectrum in the lower half plane (see 
Figure [3]). The dominated convergence theorem allows us to insert a factor L S /(L + iu) 3 , 



L 3 



— lim Im r 

7T L^oo J c (L + iuj) 



S UL ,u R (u) du 



The function L 3 / (L + iuj) 3 has no poles in the lower half plane and decays cubically for 
large \uj\. Furthermore, according to (|4.9p . the resolvent grows at most linearly for large 
\u\. This allows us to deform the contour in such a way that C is closed in the lower half 
plane on larger and larger circles \uj\ = R. In the limit R — > oo the contribution along the 
circle tends to zero. Thus we end up with the integral along the contour C ULjUR . ■ 



Our next goal is to get rid of the "convergence generating factor" L 3 / (L + iu>) 3 in (|7.2p . 
We shall use the fact that when we take the difference S ULjUR — Soo and evaluate it with 
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a test function, the resulting expression has much better decay properties at infinity (see 
Lemma HT3|) . We choose a contour Cx> which coincides with C UL)UR if Rew and 
always stays inside Q (see Figure [2]). 

Lemma 7.2 For every ^ E Cq°((ul,ur) x S 2 ) 2 , 

<*,# c *> = -1ml I <^ , S UL , UR ^> dw - f <*, 1 S" 00 ^> dw) . (7.3) 

^ V J Cuj^,u R J Coo J 

Furthermore, 

<^,E C ^> = ^2 In with (7.4) 

neIN 

In = ( / <y,QnS UL)UR ^> dw - I <^,Q n Soo^> dw) 

\ J Cuj^.uq ** Coo J 

( / <$>,Q n S UL , UR y> duj - / Q n ^oo*> dw) .(7.5) 

The series in |7.^|) converges absolutely. 



We point out that the above integrals are merely a convenient notation and are to be given 
a rigorous meaning as follows. We formally rewrite the integrals in (|7.3p (and similarly 
in (USD) as 

/ {S UL)UR - 5oo)^> dw + i - J <*,S UL , Uil ^> dw. (7.6) 

J Coo \^C U j^,u R " Coo J 

Now the first summand is well-defined according to Lemma 14.31 In the second summand, 
the integrals combine to an integral over a bounded contour, and this is clearly well-defined 
because the contour does not intersect the spectrum of H Ul ^ Ur . 

Note that in (|7.5p we cannot combine the integrals over Cy i ^^u R and Cu^^u^ (and sim- 
ilarly over Coo and Coo) to the imaginary part of one contour integral because Q n in 
general does not commute with S ULjUr , and so the integrands in (|7,5p need not be real. 
For notational convenience, we abbreviate the second line in (17. 5p by "— ccc" (for "complex 
conjugated contours"). 

Proof of Lemma \ 7. S\ According to Corollary 14.21 the resolvent Soq(w) is holomorphic for 
winQ and grows at most linearly in Thus for all L > 0, 



fa 



Scioiuj) duj = . 



'Coo ( L + 

Combining this identity with (17. 2p . we obtain the representation 
<^>,E C ^> = - lim 

7T L^oo 

x Im I f it L ■ ^ S u L ,u R ^> dw - I , T L . 5oo^> dw 

\Jc UL , UR (L + kj) 3 R J Coo (L + iwf 
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We now rewrite the integrals according to (|7.6p . If replace the contour C UL:UR by C^, the 
integrands combine, and we obtain the expression 



— lim Im f — : 

7T L^oo Jq [L + tUJ 



lim Im / lT — A% (S UL}UR - Soo)^> du . 



The estimate (|4.1(jp allows us to apply Lebesgue's dominated converge theorem and to 
take the limit L — > oo inside the integrand. The error we made when replacing C ULtUR by 
Coo is 




~ lirri^ lm{ ( / - / ] , r | , A3 <^,5 UijUB ^> du; 

Now the contour is bounded, and since the factor S UL)UR i &> is bounded, we can again 
apply Lebesgue's dominated convergence theorem to take the limit L — > oo inside the 
integrand. This gives (|7.3p . 

Note that our contours were chosen such that the condition (|2.6p is satisfied for a 
suitable constant c > 0, and so Lemma [2J] applies. Using completeness of the (Qn)neiN 
(see Lemma |2. II (iii)), it immediately follows from (|7,3p that 

= ( / XI <^iQnSu L ,u R ^> duj - I <^,Q n 5oo^> duj \ - CCC. 

\JCu L ,u R ng]N ''^ neM / 

Again replacing the contour C UL)UR by Cqo, we obtain the expression 

— / Q n (S UL UR - S' 00 )^> cL; - ccc . 

Zttz In — 



neIN 



According to (|4.1ip . the summands decay faster than any polynomial in A n . Applying the 
angular estimates (|2.8j) and ()2.9p . we conclude that the sum over n converges absolutely, 
uniformly in u G Cqq. Thus the dominated convergence theorem allows us to commute 
summation and integration, and the series converges absolutely. It remains to consider 
the expression 

~2~- [ I / ) 2 <*' QnS UL ,u R ^> du - ccc. 

\''Cu L ,u R J Ceo J ng]N 

Now the contours are compact, and thus the absolute convergence of the n-series is uniform 
on the contour. Hence we can again apply Lebesgue's dominated convergence theorem to 
interchange the summation with the integration. ■ 



We shall now deform the contours C ULjUR and Cqo and analyze the resulting integrals. 
Our aim is to move the contours onto the real axis such that they reduce to an w-integral 
over the real line. It is a major advantage of (|7.4p that the series stands in front of the 
integrals, because this allows us to deform the contours in each summand I n separately. 
Moreover, since our contour deformations will keep the values of the integrals unchanged, 
Lemma 17.21 guarantees that the series over n will converge absolutely. Thus we may in 
what follows restrict attention to fixed n. 
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For given n, we know from Section that the function <^,Q n S ULUR ^f> is meromor- 
phic, and all poles are points of a{H UL<UR ). For the integrals over in (|7.5p . we cannot 
use abstract arguments because we have hardly any information on the spectrum of 
(we only know from Lemma [4.1l that the spectrum lies outside the set f2, (|4.6p . but it may 
be continuous and complex). But from the separation of the resolvent we know that the 
operator Q\Soo is well-defined and bounded unless the Wronskian w(c/>, <p) vanishes (see 
Proposition 15,41 and (|5.7p ). If this Wronskian were zero and Imcu < 0, this would give 
rise to a solution <j) of the reduced wave equation which decays exponentially as u — > ±oo. 
Such "unstable modes" were ruled out by Whiting [21]. We conclude that <^,Q n S 00 ^> 
is analytic in the whole lower half plane {Imw < 0}. 

Using the above analyticity properties of <^f, Q n S UL ,u R ^> and <\J r , Q n Soo^>, we are 
free to deform the contours C ULtUR and in any compact set, provided that C UL:UR never 
intersects a n {H ULiUR ). In particular, choosing w_ and uj + real and outside of a(H UL;UR ), 
we may deform the contours as shown in Figure HI We let E[ ul _ iU)+ \ be the projector on all 



UJ- III 


uj + &{H ULjUr 


II 

_____ N 

IV 


\ \zj \y \y 


<y \y \y v w w \y \zj / 

Cu L ,u R J 


11 on 

1 — ^^Jv 




Coo 



Figure 4: Contour deformation onto the real axis 



invariant subspaces of H Ul>Ur corresponding to real u in the range u;_ < uj < u> + , 

QnE[u_ >ClJ+ ] = _____ Qn{^) E w . 



The next lemma shows that the integral over Cm U Cm equals Q n E\ UJ _ jUJ+ \, whereas 
the integrals over the contours II and IV can be made arbitrarily small by choosing \u>±\ 
sufficiently large. 

Lemma 7.3 For every ^ E Cq°((ul,uji) x S 2 ) 2 , n G N, and e > there are lU-,lu+ € 
M\o~(H UL ^ UR ) such that 



I n + <y,Q n E [uJ _, w+] y> 



1 

2vri 



d JC Z 



< e, 



where I n are the integrals |7.5[ j and Cj is any contour in the lower half plane which joins 
lu- with u) + (see Figure^. 

Proof. Lemma 14.31 yields that by choosing cj + and — cj_ sufficiently large, we can make 
the contribution of the contour IV arbitrarily small. The integrals over Cm and Cm 
combine to contour integrals around the spectral points on the real axis, 



1 

"2vri 



'Cm JCju 



<y,Q n S u , UR ^> duo 



-y. 



<y,Q n S UL:UR iS>> duj, 



/ J\(jJ—U}'\=§ 
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where the sum runs over all u>' £ a(H UL UR ) n [cj_,u; + ], and 5 must be chosen so small 
that each contour contains only one point of the spectrum. If we let 5 —* and use 
that Q n depends smoothly on to, one sees that the integrals over the circles converge to 
-2iri <\t, Q n E w i^f>. We conclude that 



1 

"2vri 



Cm J Cm 



It remains to show that by choosing \<*)±\ sufficiently large, we can make the integral 
over the contour arbitrarily small. According to Lemma l2.lt for sufficiently large \oj±\ 
the angular operator A w is diagonalizable for all u on the contour II. Thus we can assume 
that the nilpotent matrices M in the Jordan decomposition (|5.9p all vanish. Hence we can 
separate the resolvents according to Proposition 15.41 to obtain 

Qn {S UL ,u R ~ 5oo)^> = Yl <*' ® x AT ^> ' 

AeA„ 

where AT\ is the operator with integral kernel 

AT x {u;u',$') = (r 2 + a 2 ) 2 (s Ul>Ur - Soo ) (it, u') [ H \ ) ; , \) ; , 

K \Lop{u,v) u) a{u ,v) 

Since the functions p and a are smooth and the angular operators Q\ are bounded (|2.9p . 
it suffices to show that for every e > and g E C^((ul,ur)), we can choose w± such that 
for all w on the contour //, 

00 /»oo 



00 J —00 



du du' g{u) g(u') (s [uL)Ur] - s^){u,u ! ) < e 



Let us derive a convenient formula for S[ UL)Ufl ] — Soo- We let 0i and 02 be the two 
fundamental solutions which satisfy the Dirichlet boundary conditions 0i(ul) = = 
4>2{[ u r)- Likewise, we let and be the two fundamental solutions in infinite volume 
as constructed in Corollary 16.41 Furthermore, assume that ul < u < u' < ur. Then, 
according to (|5. 



1 1 ' v 

SU L ,U R (U,U') = 01 (u) <p 2 (u') , Sooi^u') = 4>(u) (j)(u') 

W{(j)l,<p2) W(<P,4>) 

Expressing 0i as a linear combination of and 2 , 

<f>l(u) = <j>(u) 2 (u L ) - (j)(u L ) 4>2{u) , 

and substituting into the above formula for s Ul ^ Ur , we obtain 

s UL ,u R {u, v!) = — — 1 , U(u) M^l) - 4>(u L ) M u )) <M«') 

0(u) 2 (« ) 



W(0,0 2 ) 02^l) w(0,0 2 ) 

In the first summand, we can express 02 in terms of and 0, 

02 M = 0(«) <K«r) - 0(«k) 0(") • (7.7) 
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This gives 



1 



w(<j),(j)2) 



4>(u) <h(u') 



1 



-4>{u R ) w(4>,(p) 



We conclude that 



?U L ,U R 



Soo)(u,u) 



(j){u R ) 4>(u) <p(u') , ,. 

'77 \ TOT + s ^ u ) • 

4>{ u l) <p2(u) <p 2 {u') _ 4>(u R ) (f>(u) <fr(v/) 



and because of its symmetry in u and it', this identity is also valid in the case u' < u. 
Using (|7.7f) and the notation (I6,26p . we get 



H u l) <h{v) H u r) n«) 4>W) 



</>[u L ,u R ] H u b) 0) 4>{u R ) w(<j), t 



(7.8) 



We choose |cj±| such that 



fUR 

/ ImVV(w±) € 



2Z + 1 



7T . 



According to the estimate (j6.7|) in Lemma [6.21 we can arrange that the function f™ R Im y/V 
is nearly constant on the contour II, and thus 



sin 



2 / ImVV 



> 



(7.9) 



Propositions 16. 5[ 16.61 and 16.71 allow us to estimate each term in (|7.8p by the corresponding 
term in the WKB approximation. According to (|7.9p . the factor | sin(. . .)| _1 which appears 
in Proposition 16.71 is bounded. Choosing e sufficiently small, We thus obtain the estimate 



< 2 



,)(«,«' )| 
d(tii) a 2 (u) a 2 (ii') 



exp ^2 y R ReW^ + 2 



a(u R ) d{u) d{u') 



a [u L ,u R ] a{u R ) w(d,a) 
where we introduced the function 

ot2(u) = \d(u) a(u R )\ + \d(u R ) a(u)\ . 

Using the explicit formulas (|6.24l I6.28P together with (|7.9|) , we get 

1 



d(u R ) w(d, a) 



,(7.10) 



\w(d,d)\ > \\/V(u)\ \d(u) a(u)\ , \®[u L ,u R ]\ > -^\a(u R ) d(u R )\ . (7.11) 

Substituting these bounds into (|7.10p . we get an estimate for \s UL:UR — Soo\ in terms of 
expressions of the form 

1 



exp 



2 / R ReW 



"L 



a{u\) 



d(vi) 



a[u k ) 



d(v k ) 



(7.12) 
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with Ui,Vi S [ul,ur]- The quotients of the WKB wave functions have according to (|6,10p 
the explicit form 



d(u) 



d(v) 



V(u) 



V(v) 



exp ^2 J R ReVV^j 



(7.13) 



and similarly for a. The inequality (|6.6p shows that the exponentials in (|7.12p and (|7.13p 
are bounded uniformly in u. Furthermore, it is obvious from (15, 5h that V(u)/V(v) is close 
to one if | a; | is large. We conclude that on the contour II, 

\(Su L ,u R ~ SooXm,?/)! < 



VWW\ J 

and this can be made arbitrarily small by choosing \uj±\ sufficiently large. 



We are now in the position to prove our main theorem. 
Proof of Theorem According to Lemma 17.31 

-— lim ( I ~ I ) <y,QnSooV> doj = I n + <^,Q n E R ^> , 

llll u;±-^±oo \J Ci JcjJ 

where £r denotes the projector onto the invariant subspace corresponding to the real 
spectrum of H UlUr . Here the uj± are to be chosen as in Lemma 17.31 and Cj is again any 
contour which joins uj- with to + in the lower half plane. Suppose that the contour C e 
intersects the lines Re a; = oj± in the points cj + — i8+ and lo- — i8-, respectively. Then we 
choose the contour Cj as follows, 

Ci = (u^-i[0,5-})U(C e n{u-,uj + ) + iR)\j(uj + -i[0,5 + ]). 

The first and last parts of the contour have lengths <5_ and <5+, respectively, and these 
lengths clearly tend to zero as uj± — > ±oo. Furthermore, it is obvious from Propositions 16.51 
and !6.6l as well as (|7.1ip and (|7.13p that the integrand is uniformly bounded on these parts 
of the contour. Hence the contribution of these contours tends to zero as uj± — ► dboo. Thus 



1 

~2wi 



c e J a 



QnSoo^> duj = I n + <^,Q n E R ^> 



According to Lemma 17.21 and Lemma 12.11 the right side of this equation is absolutely 
summable in n and 

(4+ <*, QuEr^>>) = {Ec + E R ) ^> . 

n 

Since the spectral projectors in the Pontrjagin space Tt UL ,u R are complete, Ec + Ej& = 1. 
We conclude that 

-^E(/ Cs "/J <*,QnS»*> d» = <*,*> . 

Polarizing, we obtain for every \F 6 Co°((ri,oo) x S 2 ) 2 the simple identity 

* = ~h £ ( L - L) du, . (7.14) 



c e J a 
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The integral and sum converge in Li ' . 

If we apply the Hamiltonian to the integrand in the above formula, we obtain according 
to Proposition 15.41 

H Q n Soo^ = (H — u) QnSoo^ + u> QnSoo^ = (holomorphic terms) + u Q n Soo^ ■ 

The holomorphic terms are holomorphic in the whole neighborhood of the real axis en- 
closed by C £ and C e (see Lemma 12.11 (i)), and therefore the contour integral over them 
drops out. We conclude that applying H reduces to multiplying the integrand by a factor 
u>. Iteration shows that the dynamics of ^ is taken into account by a factor e~ lult , 

•w = -si ~4) 

Comparing this expansion with (|7.14p . one sees that the integrand in the last expansion 
is equal to the integrand in (|7.14|) if is replaced by ^{t). Since ^f(t) is smooth and by 
causality has compact support, we conclude that the integral and sum again converge in 
Lf oc . Finally, using that the contour integrals in this formula are all independent of e, we 
may take the limit e \ of each of them. ■ 
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